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Abstract 

We compute the fermionic radiative contributions to the decay H + — > V^+(*) J 4° 
in the framework of models with two Higgs doublets (2HDM), for the case of an 
on-shell and off-shell W. We show that, in the majority of the cases, current mea- 
surements of the p parameter suggest M H ± > Ma and such decays could invalidate 
current charged Higgs searches or aid detection in the region M H ± ~ My/. We find 
that the radiative corrections may approach 50% for small values of tan (3. 
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1. Introduction 



The search for the Higgs boson (0°)0 of the Standard Model (SM)@ is one of the major 
challenges for present and future colliders. In recent years there has been growing interest 
in the study of extended Higgs sectors with more than one Higgs doublet || . The simplest 
extension is the Two Higgs Doublet Model (2HDM), and such a structure is required for 
the Minimal Supersymmetric Standard Model (MSSM). Models with two (or more) Higgs 
doublets predict the existence of charged Higgs bosons, and their discovery would be 
conclusive evidence of an extended Higgs sector. In the 2HDM extension of the SM, 
from the 8 degrees of freedom initially present in the 2 Higgs doublets, only 5 remain 
after the electroweak symmetry breaking and should be manifested as physical particles: 
2 charged Higgs scalars {H ± ), 2 CP-even scalars (h° and H°) and one CP-odd scalar 
(A ). Accurate predictions for the branching ratios (BR) of these particles are required 
in order to facilitate the searches and in this paper we consider the radiative corrections 
to the decay -> A°W ( -*\ In the non-supersymmetric 2HDM (hereafter to be called 
simply 2HDM) the masses Ma and M#± may be taken as free parameters and so one may 
consider both the case of an off-shell and on-shell W. This is in contrast to the MSSM 
in which Ma and M H ± are correlated and the two body decay is never allowed. We shall 
show that current measurements of the p parameter strongly suggest M H ± > Ma for 
M H ± > 100 GeV. 

Recently it has been shown that the decay H ± — > A°W* may be dominant or even close 
to 100% in the 2HDM (Model I) over a wide range of parameter space relevant at LEP-II 
U|. This would affect current charged Higgs searches at LEP-II 0, and the Tevatron 
[0] which only assume the decays H — > tv t and cs. We therefore feel it important to 
calculate the fermionic radiative corrections to this potentially strong tree-level process. 
An additional use of the three-body decay would be the possibility of detection in the 
difficult Mff± Pd My/ region, which is considered marginal if H decays conventionally 
to two fermions. Although a thorough analysis is beyond the scope of this paper, the 
three-body decay would give rise to high multiplicity signatures of more than 4 jets, with 
a possibility of detection above the strong WW background. We note that the 2HDM 
with the popular Model II type structure cannot possess a H ± in the discovery range of 
LEP-II due to constraints from b —>■ 57 || (see also ref.JJ which derives the lower bound 
M H ± > 165 GeV), while H ± in Model I avoids such constraints and so may be light. We 
note that it is possible to have the Model II type structure and weaken the above bound 
on Mu± in a 2HDM which relaxes natural flavour conservation (NFC) JTDJ or a general 



model with N{> 3) doublets [nj. In this paper we are concerned with the 2HDM which 
imposes NFC. Limits on M#± from the Tevatron are tan f3 dependent since one requires 
a significant BR(t — > H + b) in order to obtain a visible signal. In the 2HDM with the 
Model II type structure this BR can be significant for small (< 1) or large (> 40) values 
of tan (3. For the Model I type structure it is only possible at low tan (3. 

Current mass bounds from LEP-II for the A of the MSSM force M#± > 110 GeV in 
this model, thus taking H ± out of the LEP-II discovery range |T2|. In addition, a recent 



analysis of the MSSM charged Higgs contributions to b — > 57 |R| requires M#± > 110 



GeV, a limit valid in both the MSSM and its simplest extension by adding a Higgs singlet 
superfield (NMSSM). Therefore from the point of view of charged Higgs phenomenology at 
LEP-II one may consider the 2HDM (Model I) but not more popular extended structures. 
We will present results for the case of W on-shell and off-shell for charged Higgs masses 
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of interest at LEP-II and the LHC. The paper is organized as follows. In Section 2 
we introduce our notation and the models in question. In Section 3 we evaluate the 
fermionic one loop corrections for the case of an on-shell and off-shell W, while Section 
4 displays the counterterms. In Section 5 we present our results, and Section 6 contains 
our conclusions. 



2. Notation, couplings and lowest order results 
2.1 Notation and relevant couplings 

In this paper we will use the following notation and conventions. The momentum of the 
charged Higgs boson H + is denoted by pn (ph is incoming), p w is the momentum of the 
W + gauge boson and pa the momentum of the CP-odd ^4° (pw and pa are outgoing). 
The relevant part of the lagrangian describing the interaction of the W ± with H ± and 
comes from the covariant derivative and is given by: 



£ = ^-W+{H- d>)+h-c 



2s 



w 



(2.1) 



This interaction is model independent (SUSY or non-SUSY) and it depends only on 
standard parameters: electric charge (e) and Weinberg angle (s\y = sm 6w). 

As we are concerned with the fermionic one loop corrections, we will give hereafter 
the relevant couplings. In the 2HDM there exist four different ways to couple the Higgs 
fields to matter (we assume natural flavour conservation Rl4|). The two most popular 
are: Model I: The quarks and leptons couple only to one of the 2 Higgs doublet exactly 
as in the minimal standard model. Model II: To avoid the problem of Flavor Changing 
Neutral Currents (FCNC), one assumes that one of the 2 Higgs fields couples only to 
down quarks (and charged leptons) and the other one couples to up quarks (and neutral 
leptons). Model type II is the pattern found in the MSSM. 

In general, the couplings of the charged Higgs boson H , Goldstone G , CP-odd A 
and the gauge boson W ± to a pair of fermions are: 



H + ud 
A°uu -- 
Where: 



v l (! ~ 7s) vR (1 + 7s) 



(l-7s) 



G + ud=G L ud ^ r ^ + G ud 



<R I 1 + 7s) 
2 ' ~ ud 2 



^uu 75 



Y L - 

I ud — 



I ud 



Y 



Y, 



^ud 



A°dd = Y ddlb 



W„ud = —i — —'Jim 



I 1 - 



gV ud m u 


I ud — 


gV ud m d 


V2M W tan (3 


V2M W tan (3 


gV ud m u 


yfi _ 
I ud ~ 


gV ud m d taxi/3 


V2M W tan (3 


V2M W 


gm u 


Ydd = 


gm d 


2M W tan (5 


2M W tan (3 


gm u 


Ydd = 


m d g tan (3 


2M W tan (3 


2M W 


gm u V ud 


riR 

U ud ~ 


gm d V ud 


V2M W 


V2M W 



V2 

Model I 

Model II 
Model I 
Model II 



(2.2) 



(2.3) 
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V u d is the Kobayashi-Maskawa matrix element. It is worth noting that Models I and II 
are not very different for the top-bottom loop corrections at low tan f3 because the term 
m t / tan/5 will dominate and it is common to both types. 



2.2 Lowest order results 

The lowest-order Feynman diagram for the two body decay H + — > A°W + and for the 
three body decay H + — ► v4°iy* — ► A°ff are depicted in the following figure: 




Figl.a 




Figure. 1 



In the Born approximation, the decay amplitude of the charged Higgs into on-shell CP- 
odd Higgs boson A and the gauge boson W + (Fig.l.a) can be written as: 

M\H + ^W + A°) = ^T» where rft = i^-{p H +p A )„ (2.4) 

Here e is the W polarization vector. We then have the following decay width: 

T -= K s>S* W MU Xi{M ^ MlM2w) <2 ' 5> 

where X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) is the familiar two-body phase space 
function. Note that in the MSSM the two-body decay of the charged Higgs boson into 
iy + v4° is kinematically not allowed. 

Below threshold, and taking into account that the virtual W* decays into a pair of 
fermions //' (/ ^ t) (Fig.l.b) which we will take to be massless, the Dalitz plot density 
for this three-body decay H + — > y4°H^ + * — > A°ff is given by 



where 



dr off n OL 2 „ [(l_ Xl )(l_ X2 )_^] 



dx\dx2 327rs^/ ([1 — x\ — x 2 — + %] 2 + nwlw) 



KA,W — ~ TTo ) lW 



M 2 H± ' ,w M 2 H± > 

Tw is the total width of the W gauge boson and Xi = 2Ei/M H ± are the scaled energies 
of the massless fermions in the final state. We note that in the non-SUSY 2HDM null- 
searches at LEP in the e + e~ — > h°A°, h°Z channels eliminate regions in the M A , Mh plane 
|J, The excluded region does not have a simple shape, and there are still areas which 
allow Ma + < 90 GeV. Thus Ma may be taken as light as 10 GeV. This is in contrast 
to the MSSM in which one can derive individual lower limits on the masses, of M^ > 70.7 
GeV and Ma > 71.0 GeV [16]. Therefore the off-shell decay in the 2HDM can be relevant 
even for a small M#± (< 80 GeV) in range at LEP-II. 
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3. Fermionic radiative corrections. 



We have evaluated the fermionic radiative corrections to H + — > W + A° (for both the 
on-shell and off-shell W) at the one loop level. This set of corrections is Ultra- Violet 
(UV) divergent. The UV singularities are treated by dimensional regularization [|17] in 
the on-mass-shell renormalization scheme. 

The typical Feynman diagrams for the virtual corrections of order a are drawn in figure 2. 
These comprise the vertex correction (Fig.2.ai, Fig.2.ci2 ), W + -W~ self-energy (Fig. 2. 03) 
and the mixed W + -G~ self-energy (Fig.2.a4). Note that diagrams 2.a 3 and 2.a 4 are 
not to be considered if the gauge boson W is on-shell. These contributions have to 
be supplemented by the counterterm renormalizing the vertex H + A°W~ (Fig.2.ci), the 
counterterm for the off-shell W gauge boson self-energy (Fig.2.C2) and by the counterterm 
for the mixing W-G (Fig 2.C3). These Feynman diagrams are generated and computed 
using FeynArts and FeynCalc |18[ [19| packages. We also use the fortran FF-package 



in the numerical analysis. Note that in the general 2HDM, the vertices W + A°G~ 
W + G°H~ and A°H + H~ are not present, and so the mixing G + -H~, G°-A° and W + -H~ 
does not give any contribution to our process. 
The one loop amplitude Ai 1 can be written as: 

M\H+ -> W + A°) = e ;r^ (3.1) 

Using Lorentz invariance, r M can be projected as: 

(Thp^ + Twp^) (3.2) 



2sw 

Th and Tw can be cast as follow: 

-p -^vertex , -\^W + W~ < rW + G~ < ^vertex , xrW + W~ < X~rW + G~ (o o\ 

1 W — 1 iy + L W + 0L W + OL w + 0L W (6.6) 

T H = T^ rtex + T^ +w ~ + ST v ^ rtex + 6T^ +W ~ (3.4) 

Where T^ r jf x , T^h^ and r{y +G are respectively the contribution of the two vertices, 
the contribution of the self-energy of the W and the contribution of the mixed W + G~ 
self-energy; 5T^ r }f x , ST^^ and <5rj^ +G are the counterterms needed to remove the 
UV divergences contained in T^ r jf x , rjy^/ 17 and r|y +G . In what follows, we write the 
above one loop corrections explicitly. The expressions for the counterterms can be found 
in Section 4. 



3.1 Vertex with u— u— d exchange: Fig.2.ai 

The amplitude of the u-u-d quarks contribution to H + A°W + vertex is given by 



(ml - ml)Y u L d B (0, ml, ml) - 2M 2 H± {(mlY u L d + m D muY^ d )C + 

YuVwC, - 2C 00 + (-M\ + M 2 H± - p 2 w )C 12 - 2M 2 A C 22 ]} ) (3.5) 

= N cAfuu ( Y u L d B (M 2 H± ,ml, ml) + mu(m v Y^ + m D Y*)C - 

Y u L d {(M 2 A - M| ± )d + M\C 2 + 2C 00 + (M\ - M 2 H± + p 2 w )C u + 

(3M 2 A - M 2 H± + p 2 w )C l2 + 2M 2 A C 22 } ) (3.6) 
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with A , B , Ci and CV,- the Passarino-Veltman functions [21 which we define in Appendix 
A. Nq = 3 for quarks and 1 for leptons. All the Ci and have the same arguments: 
(Pw, M 2 H±1 M 2 Al ml, m 2 D , mfj) 

3.2 Vertex with d-d-u exchange: Fig.2.a2 

The amplitude of this diagram can be obtained from the above one just by making the 
following replacement: 

= r%*y[ mu «— > m D , Y^ d < — > Y^ d , Y uu -> Y dd ] 

The total contribution of vertex is: 

-pvertex -puud , -pddu 

1 W,H — 1 W,H "T 1 W,H 

3.3 VK + -iy^ self-energy: Fig.2.a 3 

The contribution of W self-energy Fig. 2. 03 evaluates to 



N c a 



u) 



T H = 2/2 ^2\i A ^ m u) + m D B o(Pwi m D, ™u) - 2B 22 (p w , m D , m 

zns w\Pw ~ m w) 

+PwBi(p 2 v ,m 2 Dl m 2 u )} 

r w W = — A — 2 ( N 2° a 172 A M m u\ + m 2 D B (p 2 w , m 2 D , m 2 v ) - 2B 22 (p 2 w , m 2 D , ml, 
q = 7ls w\Pw ~ 1V1 W) 



+^S 1 (p^,TO^m^) + 2(M| ± -Ml)(5 1 (^,m^,m^)+5 21 (p^,m^,m^)} (3.7) 
3.4 W-G mixing: Fig.2.a4 

In accordance with Lorentz invariance, the mixing self-energy W-G is proportional to 
and evaluates to 



r WG _ N c a(M H± M A ) 2 2 2 2 2 ™2i R( ' T1 2 _2 2\\ 

^ lvl w s w\Pw ~ 1V1 W) 



r% G = (3.i 



4. On— mass— shell Renormalization. 



The parameters entering the tree-level amplitude in eq. (|2.4j) are all standard model pa- 
rameters (e and Sw)- This fact will render the one loop renormalization rather simple, 
in the sense that all non-standard parameters appearing first at the one loop level (like 
tan/3), will not get renormalized. This is in contrast to the calculation in |22|] for the pro- 



cess H + — > hW + which explicitly contains the factor cos 2 {(3 — a) at tree-level. Therefore 
renormalization conditions related to the definition of tan f3 are not explicitly needed here. 
We will need, however, to renormalize the electric charge, the Weinberg angle, charged 
Higgs wave-function, CP-odd Higgs wave function and W gauge boson wave function. 
In our case the W ± gauge boson mixes with the Goldstone boson G , by virtue of the 
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Laurentz invariance of the self-energy; therefore T^ +w is proportional to pjf and so if 
the W is on-shell the mixing would have a vanishing contribution but in the off-shell case 
we have to take this mixing into account. 

In what follows we will follow to an extent the on-shell-renormalization developed 
by R. Santos et al [22[ which is the generalization to the 2HDM of the I. Aoki et al 

The crucial point in this scheme is that all 
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on-shell renormalization scheme 
fields and masses are renormalized after the diagonalization of the bare mass matrices. 
Another important point in this scheme is that the gauge fixing is written in terms of 
the renormalized parameters and fields and as a consequence it does not contain any 
counterterm. 



4.1 Vertex H + A°W + counterterm 

To obtain the renormalized vertex W~A°H + vertex we have to make the following sub- 
stitutions in eq. fl2.1|) : 

W, - Z^W, (4.1) 

H* - Zl£ H+ H± (4.2) 
A - Z l J 2 A° 
e — > Z e e = (1 + 5Z e )e 

~ K • $M*r (4.3) 
Mf -> Mf + 5M 2 Z 

Note that in the on-shell scheme, the Weinberg angle is defined as: sly = 1 — ^ttt- 

z 

Therefore the counterterm of s\y is completely fixed by the counterterm of the W and Z 
boson masses and is given by: 



5_sw lc 2 w SM 2 , 5 M 2 Z 
s w " 2sl v K M* v Ml 



(4.4) 



Setting Z 1 ' 2 = 1 + UZ, one obtains the following counterterm: 



8C = ^-W+{R- f A ){Uz w + Uz A o + Uz H±H± + 5Z e - ^) (4.5) 

In the on-mass-shell scheme the counterterms can be fixed by the following renormal- 
ization conditions: 

• On-shell condition for the charged Higgs boson /f , CP-odd A and the W and 
Z gauge Bosons. We choose to identify the physical mass with the corresponding 
parameter in the renormalized lagrangian, and require the residue of the propagator 
to have its tree-level value, i.e., 

5M 2 = ReZiM 2 ) and ^=- S k tf (4.6) 

ok 2 

where S(fc 2 ) is the bare self-energy of the H , A or W. 



the electric charge e is defined as in the minimal standard model [24, 25 
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Tadpoles are renormalized in such a way that the renormalized tadpoles vanish: 
Th + 5t h — 0, Th + Stjf = 0. These conditions guarantee that t>i i2 appearing in the 
renormalized lagrangian are located at the minimum of the one loop potential. 



Using these renormalization conditions and as is shown in [24|, the renormalization con 



stant of the electric charge and counterterm of gauge boson mass are given by: 

5Z e = --5Z 7 , + -—8Z Zl = --^-^ + —-^T- ( 4 - 7 ) 
5M^ = Y™ W (M^) and 5M 2 Z = S| Z (M|) (4.8) 

Where T^ w , Sf z E 7 - 7 are respectively the W, Z and photon self-energies depicted in 
Fig.2.6 3i 4 i5 , T index is to denote that we take only the transverse part. We stress at this 
stage that the fermionic contribution to the mixing EJ 2 ' (k 2 ) vanishes at k 2 = 0. 

4.2 Counter term for the W self— energy and the mixing W— G 

One obtains the counterterm for the W-W self-energy by substituting eqs (|4.1] , [4.3|) in 
the W lagrangian: 

5{WW) = i{sT - ^f^)(SMw + {M 2 w -p 2 w )5Z w ) + l P ^pL(5M 2 w + M 2 W 6Z W )(4.9) 

Pw Pw 

All the counterterms appearing in 5(W fl W u ) are fixed by the renormalization conditions 
fixed above eqs. ( |4.6| , |4.8|) . 

As we have mentioned above, W + boson and G + goldstone mix. To treat this mixing, 



R.Santos et al [22] have considered the mixing of G -H which they have renormalized 



in the following way: 

H± - Z X £ H+ H± + ^i / + 2 G+ G ± (4.10) 
G ± 4 / + 2 G+ G ± + Z£ H+ H± (4.11) 

At the one loop level Z\[ 2 = 1 + 1/28 Za and Z]j 2 = 8Zij where 5Z{j = 0(a). These four 
renormalization constants together with the counterterm mass of the charged Higgs bosons 
are fixed by imposing the on-shell condition (mass located at the pole of the propagator 
and residue equal to one) and the vanishing mixing both for T<Q+ H +(k 2 ) self-energy at 
k 2 = Mfj ± and £#+ G +(fc 2 ) self-energy at k 2 = 0. Note that the Goldstone boson receives 
its renormalized mass from the gauge fixing lagrangian. Before introducing this lagrangian 
the Goldstone boson is massless, and so the renormalization conditions imposed on the 
propagator of the Goldstone and its mixing with charged Higgs boson will be fixed at 
k 2 = 0. 

At the one loop level the renormalization constants 5Zh+h+ and 5Zg+g+ are given by 

dZ H+H+ (k 2 ) l dZ G+G+ (k 2 ) l 
dZ H+H+ = — \ k 2 =M 2 H± and dZ G+G+ = — | fc2=0 (4.12) 



Performing the replacement (4.1, 4.3 and |4.1l|) in the W gauge fixing term iM^d^W^G 



generated from the covariant derivative, one finds the following counterterm for the mixing 

W+-G-: 

5(W+G~) = *PwMw( 6 -^f + \{5Z W + 5Z G+G+ )) (4.13) 
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This completes the set of counterterms needed for our study. The renormalization con- 
stants of the wave function and the mass counterterms are given in the appendix B. 



4.3 Back to counter— terms form factors 



After the short discussion in section 4.2 about the on-shell renormalization we are using, 
we are now able to give the expressions of the counterterms ST^ r ^ x , SF^h^ ,^w + ° 
defined in eq. 



5Tw rtex = ~(SZ e - + \{5Z H+H + + 5Z A o + 5Z W )) 

•WW IV /\/f2 i J2 ™2\r/7 xn/r 2 J/2 ry__ 1 //„2 *jr2 



ST% W = {(M 2 H± +p 2 w - m 2 A )5Z w - 8M 2 , - M w SZ w }/(p 2 w - M&) (4.14) 

Iw) 



5T™ W = 2{5M 2 W + (M 2 W - P 2 w )5Z w }/(p 2 w - M 2 



f 2 



ST W G = -(M 2 H± - M 2 A ){^ + SZ W + 6Z G+G+ }/(p 2 w - M 2 W ) 



5. Numerical results and discussion 

In the previous section we have summarized the analytical formulae for the fermionic 
0(a) radiative correction to the decay H + — > H^+A . In this section we focus on the 
numerical analysis. We take the following experimental input for the physical parameters 

• the fine structure constant: a = f- = 1/137.03598. 

• the gauge boson masses: M z = 91.187 GeV, M w = 80.41 GeV and T w = 2.06 GeV 

• the input lepton masses: m e = 0.511MeV , = 0.1057GeV , m T = 1.784GeV 

• for the light quark masses we use the effective values which are chosen in such a 
way that the experimentally extracted hadronic part of the vacuum polarizations is 
reproduced $Z% : 



m d = 47 MeV m u = 47 MeV m s = 150 MeV 

m c = 1.55 GeV m b = 4.5 GeV 



For the top quark mass we take m t = 175 GeV. In the on-shell scheme we consider, 
sin 2 8w is given by sin 2 9w = 1 — -^r, and this expression is valid beyond tree-level. 

In the on-shell case it can be shown that the interference term 2ReM°*Ai l , found from 
squaring the one loop corrected amplitude \A4° + TV! 1 ] 2 , is equal to r^|A^°| 2 . Hence 
the one loop corrected width Tl n can be written as Tl n = (1 + IV)r° n , with Th being 
interpreted as the fractional contribution to the tree-level width. In the off-shell case, 
and taking the final state fermions to be massless, 2ReM°* Ai 1 is again equal to IV|.M | 2 , 
although IV now has a dependence on E\ and E2 and thus cannot be factorized out of 
the phase space integral. Therefore we define the fractional contribution to the tree-level 
width as 5T a ff, with: 

riff = (1 + SToffWff 
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Since IV does not contribute to the corrected matrix element it is evident that the W + G + 
mixing has a vanishing contribution and is given in Section 3.3 for completeness. 

We now briefly consider the constraints on the masses of the Higgs bosons that can 
be extracted from current precision measurements of p°, defined by: 

P pM% cos 2 6 W { } 

Here p in the denominator contains all purely SM radiative corrections, while p° = 1 in 
the absence of new physics. In the 2HDM there are extra contributions to p° |28j and 



Ref. p9| shows that —0.0017 < 5p° < 0.0027 at the 2a level. Imposing this condition 



and using the formulae in Ref. we plot in Fig. 3 the allowable values of M#± and 
Ma- We vary all Higgs masses up to 500 GeV and respect the current experimental lower 
limits for 5000 randomly chosen values. In Fig. 3 the triangles (points) disallow (allow) 
the decay H ± — > AW*. From the figure we can clearly see that for M#± > 100 GeV the 
vast majority of the allowed parameter space satisfies M^± > Ma, thus implying that the 
decay if ± -> AW& will be open for M#± of interest at the LHC and the Tevatron. For 
M H ± < 100 GeV (i.e. the LEP-II range) it is easier to find M H ± < M A . 



5.1 On— shell W gauge boson 

We now present our results for the case of the W boson being on-shell. There are three 
unknown parameters which determine the magnitude of the one loop corrected width T\ n : 
M H ± , M A and tan /?. This is in contrast to the decay if — > hW in which the mixing 
angle a and the mass of the heavier CP-even Higgs Boson (if) enter the calculation 
| ]22|j . We stress that this latter analysis only considered the top-bottom loops, while 
we include all the fermions corrections and find that the light fermion loops are not 
entirely negligible. Moreover, there can be significant interference among the various 
contributions, both destructive and constructive. We consider both Model I and Model 
II, which have effectively identical results at small tan /?, although differ at large tan /?. 

Let us discuss first the effect of a relatively light charged Higgs (M H ± < 250 GeV) and 
a very light CP-odd (Ma ~ 35 GeV) on Th- In Fig. 4 we plot Th in Model II as a function 
of Mh± for several values of tan/3. We note first that for a fixed value of tan/3, Th is 
insensitive to the variation in Ma when M#± is varied from 120 to 260 GeV. The peaks 
correspond to the opening of the decay ff + — > th. For small tan/3 and M H ± < 170 GeV 
the correction is rather small (~ 2%); when Mjj± > 180 GeV one can reach a correction 
of 10%. In the case where tan/? is large, the effect comes exclusively from the bottom 
quark mass and is around 10 %. 

In Fig. 5 we plot r# as a function of Ma, taking Mjj± = 570 GeV and 3 small values 
of tan f3. Since we are not considering large tan f3 this plot is relevant for both Model I 
and II. For Ma less than 300 GeV or heavier than 360 GeV the effect is about 5%. When 
Ma becomes close to 2mt a sharp peak appears and this corresponds to the opening of the 
channel A — > ti, the maximal effect in this case being around 50%. For Ma away from 
this threshold value (Ma ~ 330 — > 345 GeV) and for small tan (3 one can have a correction 
of about —14% — > —41% . As tan/? increases one quickly approaches a horizontal line 
at 3.3%. These effects are explained as follows: the ttb loop correction is proportional 
to Y uu and dominates the bbt loop correction at small tan/? because m< ^> mi, . Since 
Y uu is proportional to 1/ tan/? we can explain the tan/? dependence in Fig. 5. As tan/? 
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increases the contribution of the ttb loop weakens rapidly and the dominant contribution 
to the corrected width becomes that of the renormalized e and sw, giving a fixed value 
of Th ~ 3.3% which is very insensitive to tan/5 (note that the bbt loop in Model II is 
proportional to tan f3 - see below). We do not notice an obvious correlation between M H ± 
and IV; for the optimal case considered of tan (3 = 0.5 and Ma ~ 330 GeV, varying M H ± 
from 450 GeV to 800 GeV causes T H to fall from -18% to -27%. 

In Fig. 6 we plot T# in Model II as a function of tan/? for tan/? > 20. In Model I 
all the fermion loops decouple as tan/3 increases and one has « 3.3% for tan/? > 4. 
In Model II the bbt loop dominates with increasing tan (3 and for tan /3 > 20 the value of 
T# starts to differ from the corresponding value in Model I. Again one can find sizeable 
negative corrections, with the largest occurring for smaller Ma i.e. the closer Ma is to 
2rrib, the more on-shell the virtual b quarks are. 

In Fig. 7 we show graphically the relative magnitude of the sum of the heavy quark 
loops, ttb and bbt, compared to the sum of the remaining fermion loops (Tu g ht)- Since we 
plot only low values of tan (3 the ttb contribution dominates the bbt loop and so we label 
the sum of the ttb and bbt contributions as T U b- One can see that Tu g ht is of comparable 
strength to the heavy quark loops unless Ma is close to 2m t . In addition there can be 
constructive or destructive interference, which is shown in Fig. 7 by the sign of the ratio. 



5.2 Off— shell W gauge boson 

We now consider the case of the W gauge boson being off-shell. Since the decay H ± — > 
AW* is possible for a light H ± in range at LEP-II we shall present results for Mh± = 80 
GeV, which is also in the mass region considered problematic for detection channels 
which make use of the conventional decays H ± — ► ru T , cs. As is mentioned in the 
introduction, charged Higgs bosons of Model II are excluded from the LEP-II discovery 
range by precision measurements of b — ► sj. Our discussion will therefore be focussed 
on Model I. In the massless fermion final state limit, the WW self-energy is the only 
additional contribution to the one loop corrected width for the off-shell decayQ. The 
WW self-energy is the standard diagram and does not depend on tan/?. Hence all the 
tan/3 dependence is contained in the vertex contribution and in the case of Model I is 
enhanced when tan (3 is small. 

In Fig. 8 we plot the magnitude of the one loop corrections, ST a ff, as a function 
of small tan/3 for two values of Ma- We can see that for tan/3 > 2 one approaches a 
fixed value (~ 2%) for 5T Q ff - this is to be interpreted (as before) as the fermion loops 
decoupling, leaving a tan /? independent value which comes from the WW self-energy and 
from the renormalized e and Sw in the vertex contribution counterterms. For low tan /? 
the one loop corrections are pulled negative. Very large corrections of up to —90% are 
possible for exceptionally small (~ 0.1) values of tan/?, although such values are strongly 
disfavoured by measurements of Rb which require tan/? > 1.8 (95% c.l) for M H ± = 85 
GeV M. 



1 Note that for the W being off-shell, there are extra contributions coming from box diagrams which 



will be considered in ref. 30 
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6. Conclusions 



We have computed the Yukawa coupling corrections to the decay H + — > A !^ 4 " in the 
case of an on-shell and off-shell W gauge boson. We have included in our analysis both 
top-bottom contributions and light fermion contributions, the latter being non-negligible 
and may interfere destructively or constructively with the former. Restrictions on the 
possible values of the Higgs boson masses from considering the p parameter were also 
included and found to give in the majority of the cases M#± > Ma- In the on-shell case, 
we studied the sensitivity of the Yukawa corrections to tan/3, and found similar effects 
for small tan/5 in both Model I and Model II which can reach 50% for rriA ~ 2mj. For 
large tan/3, in Model I all the fermions corrections decouple and reach a constant value 
3.3% for tan/3 > 4; in Model II, the top mass effect is suppressed while the bottom mass 
effect is increased for tan/3 > 20, allowing sizeable corrections of 10% or greater. For the 
case of the W gauge boson being off-shell, the charged Higgs bosons in the LEP-II range 
and tan/3 not too small, the corrections are rather small and do not surpass 2%. 
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Appendix A: Passarino Veltman functions 



Let us recall the definitions of scalar and tensor integrals |2T| we use: 



A.l One point function: 

The one point function is defined by: 



.2 



(2TTp) 4 - d f , 1 

• w = J d q w^\ 

p is an arbitrary renormalization scale. 

A (m 2 ) =m 2 [l + A ] + 0(d - 4) (A.l) 
The UV divergences are contained in Aq which is given by 



A, = -^— d - lE + log(4vr) + log £ (A.2) 



note that in dimensional regularization Aq(0) = 0. 
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A. 2 Two point functions: 

The two points functions are defined by: 

n / 2 2 2n (2vr/i) 4 - d /■ 0^,/^ 

50,^(^,^0,77^0 = — / d qy- 2TT7 i V2 2T 

tir z J [q z — ttlq] [(q + pi) A — mf\ 

„ 1,. » 1 f 1 , r x 2 k 2 — x(k 2 — ml + m?) + m? — ie 

B = - A + AO - / cfeLo^ i °- V- l - 

2 jo m mi 

the derivative of Bq function is defined as: 



d 

B Q [X,ml,ml] = —B [p 2 ,ml,m 2 2 }\ p 2 =x 

note that A can be expressed in term of B 

Ao(m ) = m + m Bo(0, m , m ) 
using Lorentz invariance, we have: 

B v = PipBt 

B^v = P\ i $\ il B 21 + g^ v B 22 

A. 3 Three point functions: 

The three point functions are defined as: 

y~t (22 2 2 2 2\ _ ^ f m 1) ^M?^ 

^04i^v\PliPl2iP2-> m 0-> m l-> m 2) ~ ~^~~o / G 2lT7 I \2 2l77 I \2 2T 

in J [q - mgj [{q + p x ) 2 — m{\ [(q + p 2 ) - mgj 
where pf 2 = (Pi + P2) 2 - Using Lorentz invariance, and C^j, can be written as: 

C„ = PijA + p 2/t C 2 (A.4) 

Cll + P2 V P2 V C22 + {Pl^P2 v + P2f 1 Pl u )Cl2 (A.5) 

Appendix B: Renormalization constants 

Hereafter we give all the renormalization constants necessary to compute the counterterms 
defined in eqs [4.14|: 



B.l gauge bosons self— energies 

Let ie7^(VL^-^-+V / Ri : y 2 -) the general coupling of the gauge bosons to a pair of fermions 
/ and /'. The coefficient of — g^ v of the self-energy of the gauge boson is given by: 



+m f [m f iy L 2 + V R 2 ) - 2m f ,V L V R ]B (p 2 , m 2 f „ m))} (B.L 
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•V = Z,f = f',Z L = -^(T f -e f * w )Z R = %&, with T u = \ and T d = -± 

• V" = 7, / = /', 1l = 1r = ~e f 

• V = W, f = u and /' = d W L = W R = 0. 

The renormalization constant of the electric charge is given by: 

1 ldE?(^ 



5Z e - --SZyy - ~ ^ | fc2=0 

= l^{el[-l + B (0,mlml) + 2mlB' (0,mlm 2 d )\ (B.2) 
+e 2 u [~ + B o (0, m 2 u , ml) + 2m 2 u B' (0, m 2 u , m 2 u )]} 
The mass counterterms for the gauge boson W and Z are given by: 

5M 2 W = £5™V = M 2 W ) , 5Mf = £fV = Ml) (B.3) 

B.2 Wave functions renormalization 

The wave function renormalization constants for the W gauge boson can be obtained from 
the self-energy as: 

S7 _ dT^ w (k 2 ) 

w ~ dk 2 lfc2=M ^ 



47r J 2 s 2 i M w/3 + K - m 2 u ) 2 B (0,m 2 d ,m 2 u )/M^ 



l w°w 

- (m 2 d -m 2 u ) 2 + M^)B (M^,m 2 d , m 2 u )/M 2 w 

+ [(m 2 d -m 2 J 2 + m 2 u M^-2M^ + -2m 2 d M^}B' (M^,mlm 2 u )} (B.4) 

The renormalization constants of the charged Higgs, CP-odd Higgs and the Goldstone 
boson wave function are given by: 

SZ B+H+ = ^(-(Y u L f + Y u f)B (M 2 H± ,m 2 d ,m 2 u ) 

+([m 2 d + m 2 u - M 2 H± ](Y u L d 2 + Y u f) + 4m d m u Y u L d Y u R d ) B> (M 2 H± , m 2 d , m 2 u )) 

5Z A o = -N c ^-{Y 2 d [B (M 2 A ,mlm 2 d ) + M 2 A B' (M 2 A , m\ m 2 )] + 

Y 2 u [B (M 2 A ,mlm 2 u ) + M 2 A B' (M 2 A ,mlm 2 u )]} (B.5) 



5Z G+G+ = ^ C 2 {-(ml + m 2 u )B (0, m 2 d , m 2 u ) + (m 2 d - m 2 u ) 2 B' (0, m 2 d , m 2 u )} 
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Figure Captions 



Fig. 1 Feynman diagrams for the Born approximation to the decay H + — > A°W + *, W 
on-shell (La), W off-shell (l.b). 

Fig. 2 Feynman diagrams for the one loop corrections to the decay H + — > A°W + *: vertex 
(2.ai and 2.ai), WW self-energy (2.03), W-G mixing (2.a 4 ). Charged Higgs boson, 
Goldstone boson and CP-odd self-energies (2.61,2,3). (2. 64,5,6) WW, ZZ and 77 self- 
energies. (2.ci) is the vertex counterterm, (2.C2) W self-energy counterterm and 
(2.C5) is the W-G mixing counterterm. 

Fig. 3 Scatter plot of values of M H ± and Ma consistent with measurements of p°. Triangles 
disallow the decay -> AW^> . 



Fig. 4 T H as a function of M H ± (Model II) for M A = 35 GeV, tan/5 = 0.5, 1.0, 4 and 70. 
Fig. 5 T H as a function of M4 (Models I and II) for tan/3 = 0.6, 1.0 and 2.5. 
Fig. 6 T H as a function of tan/3 (Model II) for tan/3 > 20. 

Fig. 7 T ttb /T ng ht as a function of small tan/3 for several values of M A (Models I and II). 
Fig. 8 5T ff as a function of small tan/3 and for M A = 15, 40 GeV (Model I). 
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